Electron spin phase relaxation of phosphorus donors in nuclear spin enriched silicon 
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We report a pulsed EPR study of the phase relaxation of electron spins bound to phosphorus 
donors in isotopically purified ^^Si and natural abundance Si ("''*Si) single crystals measured at 
8 K. The two-pulse echo decay curves for both samples show quadratic dependence on time, and 
the electron phase relaxation time Tm for ^®Si is about an order of magnitude shorter than that 
for "''^Si. The orientation dependence of Tm demonstrates that the phase relaxation is caused 
by spectral diffusion due to flip-flops of the host nuclear spins. The electron spin echo envelope 
modulation effects in ^'"'Si are analyzed in the frequency domain. 
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Group-V impurities in silicon have been studied ex- 
tensively in semiconductor physics. Experimental tech- 
niques such as infrared absorption, photoluminescence, 
and electron paramagnetic resonance (EPR) have re- 
vealed detailed properties of the impurity centers. EPR 
is particularly convenient for the identification of defect 
structures since the hyperfine (hf) interaction is a sensi- 
tive probe of the spatial distribution of the electron wave- 
function. For instance, Feher and later Hale and Mieher 
applied an electron nuclear double resonance (ENDOR) 
technique to this system, and measured hf interactions 
between the donor electron spins and their neighbor- 
ing host nuclear spinJ^. These experimental works, 
together with theoretical investigations'^, have deepened 
our understanding of shallow donor impurities. 

Recently, Kane and others gave a new perspective to 
the donors in Si, a playground for solid-state quantum 
information processing, since nuclear and electron spins 
in semiconductors can be regarded as well-isolated two- 
level systems: qubits^i^. If the donor electrons are 
qubits, ^^Si nuclei that have spin-1/2 and occupy 4.67% 
of the lattice sites in natural Si ("'^'Si) are decoherence 
sources as their flip-flops produce fluctuations of the lo- 
cal fields. Indeed, ^^Si-depleted, isotopically controlled 
^*Si:P exhibited the coherence time two orders of mag- 
nitude longer than "^*Si:Bi^, demonstrating that such 
nuclear spin diluted Si would be indispensable for build- 
ing a practical Kane- type quantum computer. On the 
other hand, a study of the decoherence caused by the 
spectral diffusion arising from nuclear flip-flops requires 
a material of the opposite class, nuclear spin enriched 
Si. This novel material is also interesting because of its 
similarity to III-V materials in that the electrons are lo- 
calized in a sea of nuclear spins, and more preferable for 
our purpose owing to the negligibly-small spin-orbit in- 
teraction in bulk Si, which could otherwise contribute to 
decoherence. 

In this Rapid Communication, we report the phase re- 
laxation time Tm for P donor electron spins in isotopi- 



cally purified ^^Si and "^*Si measured at 8 K. The tem- 
perature was chosen so that Tm would not be affected by 
the spin-flip time Ti. The ground-state electron can be 
excited by absorbing a phonon if the phonon energy is 
comparable to the transition energy from the Ai ground 
state to the E or T2 excited states. When returning to 
the ground state, the electron is subject to a spin-flip at a 
certain probability. This Ti process, known as an Orbach 
process, also limits Tm over the temperature range from 
10 K to 20 KSiS. While Ti is dominated by the Orbach 
process down to 6 K, and extends exponentially with 
coolingi2iii , Tm starts to deviate from Ti and becomes 
insensitive to the temperature below about 10 K2j^. Since 
our spin echo experiments require each pulse sequence to 
be repeated at time intervals much longer than Ti, we 
found 8 K to be an appropriate temperature, low enough 
for Tm not to be limited by Ti but high enough to ensure 
a reasonable measuring time. 

A Cz-grown single crystal of ^^Si, enriched to 99.23%, 
had a rectangular shape with its long axis in the [110] 
orientation. The sample contained 1.8 x 10^^ P/cm'^ 
with the compensation of 1.0 x 10^^ B/cm'^. Further 
information on this crystal is provided in Ref. IT^ A 
natgj sample was cleaved from a commercial high-quality 
wafer containing 0.8 x 10^^ V/ctc? with a negligible 
amount of compensation. The net donor concentra- 
tions of both samples were kept low so that the dipo- 
lar or exchange interactions between donors would be 
suppressecli2ii4. Pulsed experiments were carried out us- 
ing a Bruker Elexsys E580 spectrometer, and samples 
were kept in an Oxford ER4118CF cryostat. Tempera- 
ture was controlled with an Oxford ITC503 temperature 
controller. The echo-detected EPR spectra, in which the 
intensity of the Hahn echo was measured as a function 
of the external magnetic field, consisted of two Gaussian- 
shaped lines separated by 4.2 mT. The splitting is due 
to the hf interaction with '^^P, and each line is inhomoge- 
neously broadened by the surrounding ^^Si nuclei. The 
linewidths (FWHM) are 0.26 mT for "'^^Si and 1.2 mT for 



2 




r [(lO'^sf for ""'Si & (10"'s)^ for "'Si] 



FIG. 1: The two-pulse electron spin echo decay curves as a 
function of ai — 0° and 50°. Note that the unit of the 
horizontal axis differs for each sample, and the scale for each 
8. The dotted lines in '^^Si are the fits to the echo envelope 
decays. The oscillation observed in ^^Si is ESEEM. See text. 
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FIG. 2: The orientation dependence of Tm for "''*Si and ^^Si 
(left axis). The open squares (□) represent Tm for "^*Si, and 
the open circles (O) for ^®Si. Note that the unit of the vertical 
axis differs for each sample. The inverse square root of A/2 
calculated based on the method of moment is also shown by 
a dashed line (right axis). 



^^Si. In the following experiments, the external magnetic 
field was set to the center of the line at higher fields {Bq 
= 348 mT). Ti was measured using an inversion recov- 
ery method (7r-i-7r/2-T-7r-r-echo), and is 16 ms for "^'Si 
and 4.4 ms for ^^Si. As the temperature dependence of 
the Orbach process is given by l/Ti = i?exp(— A/fcT), 
where R is the rate constant and A the valley-orbit split- 
ting energy, the difference in Ti between samples could 
arise in part from a slight difference in the actual sample 
temperatures. The isotope shift of A is unlikely to cause 
the difference in Ti, since it was not observed in our in- 
frared photoconductivity measurement on the ^^Si crys- 
tal within the resolution used^S. The presence of com- 
pensation and dislocation (10^ cm~^) in the ^^Si crystal 
can alter Ti, changing Ri^ilL. However, we can conclude 
that this difference in Ti has little effect on the difference 
in Tm presented below, based on the previous assumption 
that Ti does not contribute to Tm- 

The phase relaxation was investigated using a two- 
pulse spin echo method (7r/2-r-7r-T-echo, where the in- 
terpulse delay r was varied in 800 ns steps for "'^^Si and 
40 ns steps for ^^Si. The 7r/2 pulse was 16 ns.). The 
samples were rotated around the [110] axis perpendicu- 
lar to Bo- We define 9 as the angle between Bq and [001]; 
therefore, 9 = 0° when Bq \\ [001], 9 = 55° when Bq \\ 
[111], and 9 = 90° when Bq \\ [110]. Since the echo- 
detected EPR spectra were independent of the crystal 
orientation and no other EPR signals were found, the 
alignment of the crystal from an EPR signal was not ap- 
plied here. We estimate the uncertainty in 9 to be less 
than 5°. Figure ^ shows the echo decay curves at 9 = 0° 
and 50°. Although so-called electron spin echo envelope 
modulation (ESEEM) obscures the echo envelope decays, 
they clearly obey a quadratic decay law, expressed as 
exp(— TOT^). A single-exponential term exp(— 26r) is, if 
present at all, quite small. Thus Tm can be defined as the 



time at which an echo envelope damps to 1 /e of its initial 
value, i.e., Tm = 2m~^^'^. We note that our temperature 
setting and assumption on the Ti effect are justified pos- 
teriori by the fact that the echo decay curves are not 
single-exponential and that Tm for each sample is much 
shorter than the respective Ti. 

The orientation dependence of Tm given in Fig. [21 
shows Tm to be longest at 6* = 0°, and shortest around 9 
= 50° ~ 60°—. The dependence manifests the fact that 
the phase relaxation is caused by ^^Si nuclei mutually 
coupled via the dipolar interactions. This can be verified 
by calculating the second moment M2 of the ^^Si nuclear 
spin system. M2 calculated with Van Vleck's method of 
moment is the sum of squared dipolar fields produced 
by the nuclei-'^'', and its inverse square root is a conve- 
nient measure of the nuclear T2. (M2)~^/^ for a 100% 
^^Si crystal is shown in Fig. |21 by a dashed line. Cor- 
relations between the electron Tm and the nuclear M2 
are apparent. As M2 directly reflects the strength of the 
nuclear dipolar couplings, its orientation dependence is 
understood qualitatively as follows: When Bq is along 
[111], one of the four nearest-neighbor bonds of the Si 
atoms is parallel to Bq, and this pair of nuclei gives rise 
to the strongest coupling; hence, M2 takes its maximum. 
With Bq along [001], all the dipolar couplings between 
nearest neighbors are frozen since the angle between Bq 
and the vector connecting the nearest neighbors is a so- 
called magic angle; hence, M2 takes its minimum. In fact 
such an experimental T2 has been reported for NMR of 
^■^C diamond, a material similar to ^^SiiS. As the line- 
shape studies of NMR spectra for ^'^C diamond and ^^Si 
have revealed that they share essentially the same line- 
broadening mechanism, T2 for ^^Si will show the same 
tendency as that for ^^^C diamond if measuredi2i2£. 

Although the comparison with M2 works qualitatively, 
it provides little information on the actual value of Tm- 
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Theoretical estimation of Tm must take the hf interac- 
tion between the electron and host nuclei into account as 
well as the nuclear dipolar coupling. Generally, to char- 
acterize a system where the electron phase relaxation is 
caused by the spectral diffusion due to flip-flops of the 
host nuclear spins, the diffusion barrier that prevents the 
flip-flops within its bounds must be considered^.. As the 
Fermi contact hf interaction, which is proportional to the 
density of the electron wavefunction |\l/(ri)|^, varies from 
site to site, a flip-flop of a certain pair of nuclei occurs 
only when the difference of the hf interaction between 
the pair is small enough to satisfy the condition of energy 
conservation. The condition must be evaluated for each 
pair, since (ri)p does not decrease monotonically with 
increasing but oscillates due to the multi- valley nature 
of Si. Such a theoretical treatment has been proposed 
by de Sousa and Das Sarma^, it is therefore interest- 
ing to compare our results with theirs^^ Theory predicts 
the observed angular dependence correctly, but overesti- 
mates Tm by about a factor of 3 for both samples. This 
already-reasonable agreement becomes even better if we 
take the ratio of Tm between the samples. Indeed, the 
theoretical ratio of Tm for "'^'Si to that for ^^Si falls be- 
tween 11.2 and II. 8, while the experimental ratio lies 
between 1 1. 2 and 14.4. Given the difficulty in determin- 
ing the precise Tm due to ESEEM, their calculation is 
in good agreement with our experiments. Another com- 
parison is to take the ratio of the maximum {0 — 0°) 
and minimum {6 ^ 55°) Tm- Theory yields 2.7 for "'''Si 
and 2.9 for ^^Si, compared with experimental values of 
2.0 for "'^'Si and 2.1 for ^'^Si. The larger values in the 
theory may indicate the presence of an isotropic contri- 
bution to Tm, but it is not clear at this stage whether 
other decoherence mechanisms must be incorporated or 
an improved theory of nuclear-induced spectral diffusion 
suffices to explain the discrepancies revealed here. Incor- 
porating non-Markovian nuclear flip-flop processes would 
certainly be an interesting refinement of the theory, while 
the stochastic treatment proved valid even for ^^Si. It is 
also noteworthy that recent multiple-pulse NMR stud- 
ies in Si provide a glimpse into the complicated behavior 
of this seemingly-simple dipolar coupled systemSiSiS^. 
Clearly, more experimental and theoretical investigation 
is necessary for a full understanding of the phenomena. 

We now turn our attention to the remarkable feature 
of the decay curves: ESEEM. The origin of ESEEM can 
be described briefly as follows: If the nuclear spin feels, 
in addition to the external magnetic field, the moder- 
ate hf field produced by the electron spin, the nuclear 
spin precesses around an effective magnetic field which 
is tilted from the external magnetic field, i.e., to/ is 
no longer a good quantum number. Due to this state 
mixing, formally forbidden nuclear-spin-flip transitions 
(Ams — ±I,Am/ — ±1) can occur, and interfere with 
allowed transitions to produce beats in the electron spin 
echo envelope. In two-pulse experiments for an S" = 1/2, 
/ = 1/2 spin system, the modulation contains the EN- 
DOR frequencies and and their sum and differ- 
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FIG. 3: (a) The frequency domain ESEEM spectra in ^^Si 
at 9 — 0° , 50°, and 90°. The vertical axis is shown in an 
arbitrary unit and shifted for clarity, (b) The orientation 
dependence of the ENDOR frequencies around 3 MHz. They 
show a {lll)-axis pattern. The solid (dashed) lines are for 
1/+ (i'-)- Which lattice site produces lines A, B, and C is 
explained in the text. 



ence Vj^.^v-. When many nuclei are coupled to the same 
electron spin, some combination frequencies are also con- 
tained since the two-pulse ESEEM is the product of in- 
dividual modulation functions. 

We analyzed the ESEEM spectra in the frequency do- 
main. Although ESEEM was also observed in "'^'Si, we 
treat only the case of ^^Si here because the larger mod- 
ulation depth in ^^Si facilitated the analysis. Also, the 
modulation depth is strongly angular dependent (Fig.Q, 
since the degree of state mixing depends on both the 
position of each nuclear spin and the orientation of the 
external magnetic field. To obtain a frequency domain 
spectrum, the slowly decaying part of a time domain 
spectrum was subtracted first, then the remaining modu- 
lation was Fourier-transformec^. Figure I^Ja) shows the 
frequency domain spectra at 9 — 0°, 50° , and 90° . Peaks 
around 3 MHz are the ENDOR lines, and their angular 
dependence is shown in Fig. I^Jb), from which we see a 
(Ill)-axis pattern^. The ENDOR frequencies for an ax- 
ially symmetric hf tensor with an isotropic g-factor are 
given by^^ 



Oiso + ^(3 cos^ — I) ^ ^ f3bsin2ipi^^ 



where vi is the nuclear Larmor frequency. 



the 



isotropic hf coupling constant, b the anisotropic hf cou- 
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pling constant, ipi the angle between Bq and the unique 
axis of the hf tensor. Vj is calculated to be 2.94 MHz as 
the gyromagnetic ratio of ^^Si nuclei is 8.46 MHz/T. v± 
calculated with Oiso = 570 kHz and 6 = 681 kHz agree well 
with the experimental results, as shown in Fig. |2Ib). In 
comparison with hf constants obtained from previous cw 
ENDOR experiments^, the observed peaks are assigned 
to shell E (111), i.e., four nearest neighbors of the donor. 
Lines A and B originate from (111) and (111) sites, re- 
spectively. Line C is doubly degenerate, since (111) and 
(III) sites locate each other at plane symmetric positions 
with respect to the (110) plane. The experimental data 
corresponding to line C at 6* = 0° and 10° split, however. 
This suggests that the sample was not exactly rotated 
around the [110] axis but slightly off-axis, most likely 
due to a small miscut of the crystal. This assumption is 
supported by the fact that the ESEEM in "'"^Si at 6* = 0° 
did not split (not shown). The strong peak at 5.9 MHz 
is the sum frequency, but signals from shell A (004) are 
overlapped. The fourth harmonic is also observed at 11.8 
MHz, and the third harmonic is barely visible around 9 
MHz. We did not observe the third and fourth harmonics 
in "^*Si. We also observed tiny peaks around 5.2 MHz 
throughout the angles tested. They are assigned to shell 
B (440), but the detailed angular dependence was un- 
traceable. A three-pulse stimulated echo method would 
be suitable for a more detailed ESEEM study. From 
the viewpoint of quantum computing, ESEEM clearly 
leads to quantum-gate errors. For this purpose, time do- 



main analysis is highly desirable as recently simulated by 
Saikin and Fedichkin-^. 

In conclusion, we have measured the phase relaxation 
time Tm of P donor electron spins for "^*Si and ^^Si at 8 
K. Tm for ^^Si is an order of magnitude shorter than that 
for "^'Si due to much more frequent flip-flops of the host 
nuclear spins. The orientation dependence of Tm agrees 
qualitatively with (M2)"^/^ for a 100% ^^Si crystal cal- 
culated with the method of moment, and quantitatively 
with the theory of de Sousa and Das Sarma. Frequency 
domain analysis revealed that ESEEM effects originate 
mainly from the hf interactions between the donor elec- 
tron and its nearest neighbor nuclei, as suggested by 
Saikin and Fedichkin. Our results also provide insights 
into the localized electrons in III-V materials, such as 
GaAs, whose lattice sites are full of nuclei with non-zero 
spin. Their phase relaxation would be severely controlled 
by nuclear-induced spectral diffusion, therefore the ex- 
perimental conditions must be arranged carefully so that 
the effects of nuclear spins may be suppressed, e.g., high 
magnetic fields, decoupling pulses, etc. In the near future 
we plan to prepare a series of samples with different ^^Si 
isotopic composition. Such samples will allow us to carry 
out systematic relaxation time studies of the electron and 
nuclear spins as a function of ^®Si concentration. 
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hence the Nyquist frequency is 12.5 MHz. As the data were domain spectra, 

taken from t = 320 ns, all the modulation components that 
decayed within 320 ns cannot be recovered in the frequency 



